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When systems that can undergo phase sepa- 
ration between two coexisting phases in the 
bulk are confined in thin film geometry between 
parallel walls, the phase behavior can be pro- 
foundly modified. These phenomena shall be 
described and exemplified by computer simula- 
tions of the Asakura-Oosawa model for colloid- 
polymer mixtures, but applications to other soft 
matter systems (e.g. confined polymer blends) 
will also be mentioned. Typically a wall will 
prefer one of the phases, and hence the composi- 
tion of the system in the direction perpendicular 
to the walls will not be homogeneous. If both 
walls are of the same kind, this effect leads to a 
distortion of the phase diagram of the system in 
thin film geometry, in comparison with the bulk, 
analogous to the phenomenon of "capillary con- 
densation" of simple fiuids in thin capillaries. In 
the case of "competing walls" , where both walls 
prefer different phases of the two phases coexist- 
ing in the bulk, a state with an interface parallel 
to the walls gets stabilized. The transition from 
the disordered phase to this "soft mode phase" 
is rounded by the finite thickness of the film and 
not a sharp phase transition. However, a sharp 
transition can occur where this interface gets 
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localized at (one of) the walls. The relation of 
this interface localization transition to wetting 
phenomena is discussed. Finally, an outlook to 
related phenomena is given, such as the effects 
of confinement in cylindrical pores on the phase 
behavior, and more complicated ordering phe- 
nomena (lamellar mesophases of block copoly- 
mers or nematic phases of liquid crystals under 
confinement). 



1 Introduction 

The current interest in the construction of 
nanoscopic devices [H [21 [Sj HI [5] demands a 
better understanding of the phase behavior of 
fiuids confined in pores or slits of nanoscopic 
linear dimensions P El El [TDl El [12] ■ Knowl- 
edge on the phase behavior of confined fiuids 
is a prerequisite to understand their dynamics 
[T3l [T4l ITSj . as well as for the analysis of fiow 
through very thin capillaries [HI [IT] , nanoscale 
capillary imbibition [THl [19], and related mi- 
crofiuidic or nanofiuidic devices. 
Obviously, an interplay must be expected be- 
tween surface effects on the fluid due to the con- 
flning walls, such as adsorption [201 [211 [221 [23], 
formation of wetting (or drying) layers [2H [23 
[23 [2Z], and flnite size effects [211 [211 [30] due to 
the flnite width of the capillary. However, un- 
derstanding the nanoscopic conflnement of real 
fluids consisting of small molecules is very diffi- 
cult due to additional effects, resulting from the 
lateral variation of the wall potential caused by 
wall roughness or even the atomistic corrugation 
[an [32] of the wall. 
While there has been an enormous activity to 
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study theoretically and by computer simula- 
tion confinement effects on simple fluid mod- 
els such as the Ising lattice gas model [llj or 
simple Lennard- Jones systems [lOl [12] and rich 
predictions from phenomenological theories are 
available as well [3 El El [IHl [IB [H E31 El 
[351 [361 [33 [381139], it is difficult to find perti- 
nent experiments to which such work could be 
compared. However, it is much more promis- 
ing to study confinement effects on soft matter 
systems: due to the mesoscopic length scales 
of the particles that one encounters when one 
studies mixtures of polymers and colloids [ID] 
or polymer blends [HI [12], effects due to the 
atomistic corrugation of the walls are much less 
important; also the large size of colloidal parti- 
cles enables more detailed experimental obser- 
vations; e.g. individual particles can be tracked 
though real space in real time using confocal mi- 
croscopy [l3] and interface fluctuations in mix- 
tures of colloids and polymers can be directly 
observed |111[1S]. Moreover, colloids are model 
systems for the study of phase behavior, since 
by changing suitable parameters the strength 
and range of effective interactions can be varied 
over a wide range [IHl [13 HH]- Also the inter- 
action of colloidal particles with the conflning 
walls can be tuned, e.g. by coating the wall with 
a polymer brush [l9l [501 Ell E2J and controlling 
the polymer- wall interaction via variation of the 
grafting density and/or chain length of the an- 
choring flexible polymer [S21 ES]- In particular, 
for colloid-polymer mixtures both the radius of 
the (spherical) colloidal particles and the size ra- 
tio between colloids and polymer coils controls 
the location of the critical point where the phase 
separation in a colloid-rich and a polymer-rich 
phase sets in [40] . 

Similarly, also blends of long flexible polymers 
are a very suitable model system to study the 
effect of conflnement in a thin fllm geome- 
try on phase separation experimentally as well 
[m [121 [Ml- Again, already the location of the 
critical temperature of phase separation in the 
bulk can be varied over a wide range, by suit- 
able choices of the polymeric species, and their 
chain lengths [551 [SHI [57] . In addition, charac- 
teristic lengths of the problem such as the cor- 
relation length of composition fluctuations [2Z] , 



the (intrinsic) interfacial width between coex- 
isting phases [57], etc., are much larger than 
interatomic distances, and hence also for these 
systems experimental probes are available which 
would lack sensitivity for small molecule sys- 
tems. E.g., for the study of the anomalous 
broadening of interfaces depending on the fllm 
thickness ^T\ 1^ . which is one of the char- 
acteristic signatures of the "soft mode phase" 
[211 [37] in a system with "competing walls" [55] . 
a nuclear-reaction based depth proflling method 
[51] was used. This method can resolve the 
very wide interfaces in such soft matter systems, 
while it would be unsuitable for the much nar- 
rower interfaces in mixtures of small molecules. 
As is evident from this introductory discus- 
sion, and the extensive literature that has al- 
ready been quoted, the subject is extremely 
rich, and comprehensive coverage could flU a 
whole book. Therefore the scope of the present 
review necessarily must be more narrow. We 
shall focus in this review almost exclusively 
on conflnement effects of colloid-polymer mix- 
tures [Sa [53 [601 [611 [621 [63]. Only spher- 
ical colloidal particles shall be discussed, al- 
though related phenomena can be studied also 
for mixtures of polymers with rod-like colloids 
[54] . Although extensive work has been done for 
models of polymer blends, both using the self- 
consistent fleld theory [651 [SSI [SZ] and simula- 
tions [SSI [Sa [SSI [To] , we shall not consider this 
work here, but draw attention to recent reviews 
[TT] [72] . Also, we shall not attempt to review 
the theory of wetting phenomena [211 [2S1 [2S1 [27] 
and scaling theories of capillary condensation 
[TTl [33] and interface localization transitions 
[SSI [S3 [SSI [70] but rather refer the reader to 
another thorough review [TT]. 
In Sec. 2, we shall briefly recall work on phase 
behavior of the model of Asakura and Oosawa 
(AO) [73j and Vrij [71], where colloids simply 
are described as hard spheres which may nei- 
ther overlap with each other nor overlap with 
polymers, while the latter may overlap with 
each other with no energy cost, in the bulk 
[7S1[7S1[73[781[791EQ1EI1E21ES]. In Sec. 3 
we shall discuss the phase behavior of this AO 
model when it is conflned [5S1 [SH [501 [SB [S2] 
between symmetrical walls a distance D apart. 



paying attention to the shift of the critical point 
as a function of film thickness D, and to the 
change of the critical behavior. Sec. 4 then de- 
scribes the behavior encountered for asymmet- 
ric walls [53], where it is also shown that by 
variation of the conditions at the walls one can 
gradually crossover from this interface localiza- 
tion transition to a transition which is of capil- 
lary condensation type. Sec. 5 then presents 
a summary of the results reviewed here, and 
gives an outlook on related findings in other 
systems, as well as to more complicated phe- 
nomena where the order parameter characteriz- 
ing the transition is not a simple scalar quantity 
(as it is for gas-liquid or liquid-liquid type phase 
separation) . 

2 Liquid-liquid demixing 
for the Asakura-Oosawa 
(AO) model in the bulk 

In the AO model [731 El] colloids are described 
as hard spheres of radius Re-, and hence the 
potential between two colloidal particles at dis- 
tance r from each other is 



f/n, 



oo 



(r <2Rr), Ur, 



(else) 



Similarly, polymers are described as soft spheres 
of radius i?p. Remembering that long polymer 
chains with A^ subunits have a radius i?p oc A^*^ 
with V 
Rp oc N^^"^ in Theta solvents 



0.59 in good solvent conditions [56] or 

I], the density 
Pp = N/Rp of monomers of a chain inside its 
own volume is very small, and hence polymer 
coils can interpenetrate each other with a free 
energy cost of a few fc^T (with A;^ the Boltz- 
mann constant and T the temperature) [84] • In 
the AO model, this free energy cost is neglected, 
and the polymers are treated like particles in an 
ideal gas, f/pp(r) = irrespective of distance. 
But, of course, polymers cannot penetrate into 
the colloidal particles, and hence 

[/pc(r) = oo{r<R^ + Rp), Upc{r) = (else). 

(2) 
As is well-known [10], the polymers cause an 
(entropic) depletion attraction between the col- 



loidal particles, and as a result, an entropy- 
driven phase separation occurs, if the volume 
fractions r]c, rjp of colloids and polymers are suf- 
ficiently high (Fig. [1]). Here rjc, r/p are defined 
in terms of the volume V of the system and the 
numbers of colloids and polymer, N^ and A^p, 
respectively, by 



Ajr 

r/, = —RlNjV 
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RlNp/V, (3) 



and Re = 1 will henceforth be chosen as unit 
of length. Since both rjc, rjp are densities of 
extensive thermodynamic variables, it is useful 
to carry out a Legendre transform to an inten- 
sive thermodynamic variable, where the chemi- 
cal potential fip of the polymers [or their fugac- 
ity Zp = exp(/ip/A;BT)] is used. It is customary 
to use instead of /ip or Zp the so-called "polymer 
reservoir packing fraction" rj^, 
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(4) 



Eq. (jlj) would be just the volume fraction of 
polymers in the absence of any colloids, since 
such a system simply is an ideal gas of polymers. 
It is clear that the model defined by Eqs. ([1]), ([2]) 
is a drastic simplification of reality, but in qual- 
itative respects it is remarkably accurate |4Uj . 
While various more realistic extensions of the 
AO model have been considered [79l HH [851 ESI 
EH EH Ea Eni EH [92], and sometimes better 
agreement with experiments [lll[nS] is obtained, 
we disregard such extensions here because in 
practice there are many additional effects (such 
as charges on the colloidal particles [94^ ^5\ , ad- 
sorption of polymers on the colloids [96], etc.) 
that make a quantitative comparison with ex- 
periment elusive. 

In early simulation work on the AO model 
[771 [78] a wider range of volume fractions ric, Up 
(and a much wider range of rjp) was studied, but 
only a much more limited accuracy than shown 
in Fig. [Tjwas obtained. On the basis of this work 
[771 [TSj . it was concluded that the agreement 
between simulations and the mean-field theory 
of Lekkerkerker et al. [75] is excellent. Fig. [1] 
demonstrates, however, that the relative devia- 
tion between the actual value for 77!, at the criti- 
cal point, r7pj,j. = 0.766±0.002, deviates from its 



mean-field prediction [7S] by about 30%. This 
deviation, in fact, is relatively larger than corre- 
sponding deviations between mean field theory 
and accurate simulation results for lattice gas 
models [HZ], Lennard- Jones fluids [SB], etc. In 
retrospect, this large deviation between mean- 
field theory [75] and accurate simulation results 
for colloid-polymer mixtures [801 Ell 112] is not 
surprising, since on the length scale of a colloidal 
particle the depletion attraction has a very short 
range, and the large absolute size of colloidal 
particles in this context is not relevant: it would 
be wrong to infer that colloids should behave 
mean-field like. 

Being interested in the changes in phase behav- 
ior due to confinement between walls that are 
a distance D ^ Rc{= 1) apart, relatively small 
changes must be expected, of course. For an 
analysis of these changes, and in particular for 
a study how bulk behavior in the limit D — > oo 
is approached, a very good accuracy of the sim- 
ulation data is absolutely crucial. Thus, it is 
worthwhile to briefly recall how results such as 
those shown in Fig. [1] can be obtained, since the 
methods for the study of the confined systems 
[6T1 [621 [63] are closely related to those used in 
the bulk [Hni[SIl[H2]. 

We start this recollection by emphasizing that 
for studying liquid-vapor type phase equilibria 
the grand-canonical {fiVT) ensemble of statis- 
tical mechanics is the best choice [98], since 
it avoids problems due to slow relaxation of 
liquid-vapor interfaces that hamper the use of 
the canonical ensemble [991 HOOj . Also near 
the critical point the problem of critical slow- 
ing down [lOlj is somewhat less severe in the 
grand-canonical ensemble [981 EH HOOj . and 
the inevitable finite size effects are relatively 
easy to handle by finite size scaling methods 
[211 [291 [301 [M [TO3], unhke the popular Gibbs 
ensemble [1041 1105] . So the task of the simula- 
tion is to vary the chemical potential fi of the 
colloids at fixed rj^ (as the phase diagram. Fig. [H 
suggests, 7]^ is analogous to inverse temperature 
for ordinary vapor-liquid type transitions [98] . 
where vapor-liquid phase separation is driven by 
enthalpic rather than entropic forces. Of course, 
in thermal equilibrium the average colloid frac- 
tion (rjc), which is the variable thermodynam- 



ically conjugate to fi (apart from a normaliza- 
tion factor, see Eq. [3]), increases monotonously 
with /i even when the two-phase coexistence 
region is crossed, and in the (rjc) vs. /i curve 
hence no singularity shows up for any finite lin- 
ear dimension L: only in the thermodynamic 
limit (where L —>■ cxd) this "isotherm" devel- 
ops at /i = ficocx a perpendicular part, where 
(rjc) jumps discontinuously from rj^ (vapor) to 
r]^ (liquid). However, nevertheless phase coex- 
istence is easily recognizable also in a finite vol- 
ume simulation, when the colloid volume frac- 
tion distribution P{ric) is sampled [98 ] [99 ] [TOO] . 
In the regime r]^ < {rf) < rj^, P{r]c) has a dou- 
ble peak structure, and for /i = /icoex both peaks 
have equal weight ("equal area rule" |106[ll07j ). 
In order to carry out this program, two obstacles 
need to be overcome: (i) in order to sample the 
relative weights of the two peaks of P{ric), the 
peak near rj^ representing the vapor-like phase 
of the colloid-polymer mixture and the peak 
near r^^, the liquid-like phase, the system needs 
to cross many times a region of very low prob- 
ability near rj^ = (r/^ + Vc)/'^- This problem, 
however, can be very efficiently solved by suc- 
cessive umbrella sampling |108j . Fig. [2] shows, as 
a typical example, distributions -P(?7c) that span 
almost 30 decades, (ii) The second obstacle is 
the fact that the polymer volume fraction, in the 
polymer-rich phase, can be very high (exceeding 
unity, since the polymers are allowed to overlap 
with no energy cost). Insertion of a colloid par- 
ticle at a randomly chosen position, which is one 
of the Monte Carlo (MC) moves that one needs 
to carry out in grand-canonical Monte Carlo 
simulations, almost always will be rejected: so a 
naive implementation of a grand-canonical MC 
simulation for unfavorable parameters is bound 
to fail utterly. However, this problem also could 
be overcome, by the invention of a composite 
MC move, where in a spherical region with some 
properly chosen radius r^ a randomly selected 
chosen number n^ of polymers is taken out and 
only then insertion of a colloid is attempted (the 
reverse move also exists and is constructed such 
that the detailed balance principle [991 HOOj is 
fulfilled) [SI1[S2]. 

We now return to the observation of Fig. [21 that 
a high free energy barrier AF (choosing units 



where kBT = 1) exists, which is independent of 
r)c in a broad regime of rjc around the compo- 
sition of the rectihnear diameter rjd. The inter- 
pretation of this fact is that the system in this 
region is in a state with two domains, separated 
by two domain walls, oriented perpendicular to 
the ^-direction, and connected into itself by the 
periodic boundary conditions. This is also con- 
firmed by direct inspection of the configurations 
of the system. Hence |109j 



AF = 2L27LV, L 



CX3 



(5) 



where 7lv is the interfacial tension between 
liquid- and gas-like phases, and L^ the inter- 
facial area. Thus, estimating AF for a se- 
ries of cross-sectional areas L^ of the simula- 
tion box and extrapolating the result for 7lv 
to the thermodynamic limit has become a stan- 
dard method for the MC estimation of interfa- 
cial free energies [99l llOOj . Fig. [3] shows typ- 
ical results for the reduced interfacial tension 



plotted vs. the order parameter rjl 



V^ 



and 



compares them to density functional theory pre- 
dictions [110] . These simulation results for 7lv 
are also consistent with a capillary wave anal- 
ysis [83]. Note that the coexistence densities 
r]^ , 1]^ do approach the predictions from mean 
field theory rather fast (Fig. [T]), for r^p > 1.0 the 
differences are practically invisible, however, no 
such convergence is seen for the interfacial ten- 
sion (Fig. [3]). The reason for the strong discrep- 
ancies in Fig. [3] is not clear. 
We now comment on the treatment of finite size 
effects. If one naively would take the values of 
rjc where P{ric) has its two peaks as estimates for 
7]^ and 7]^ also in the critical region, one obtains 
results as shown in Fig. HI For r]^ > 0.79 these 
estimates are independent of the linear dimen- 
sion of the simulation box, but for rj^ < 0.79 sys- 
tematic finite size effects appear. E.g., for t]^ = 
0.76 the difference r]^ — rf^ decreases systemati- 
cally with increasing L. While for 7]^ > rj^ ^^ this 
difference for L —>■ oo converges to a nonzero 
result, for rj^ < r^p,,^. it ultimately vanishes. 
While a naive inspection of Fig. Hjdoes not allow 
to estimate ^pc^ such an estimate can be ob- 
tained reliably from finite size scaling methods 

[SHiEEEniEziEHiEsiEniiinniiMiing. choos- 
ing fi = /icoex(^p) from the equal area rule, as is 



done in Figs. [H IH we define an order parameter 
m as m = T]c — {t]c) and define moments (|?ti|'^) 
{k being integer) from the distribution P{ric), 



(|m|'=) = J d7]c \m\^P{T]c). 



(6) 



Defining then the fourth order cumulant t/4 as 

[grifTus]. 

f/4 = {mY/{m') (7) 

we can invoke the result that f/4 tends towards 
unity for rj^ > tjI^^^ as L ^ 00, while U^ tends 
to 1/3 for r/p < rf^^r^ since ultimately the distri- 
bution P{i]c) in the one-phase region must be- 
come a single Gaussian centered at (rjc) [103] . 
For r/p = r/pcr, however, f/4 tends to a nontrivial 
but universal value f/4 (f/4 ~ 0.629 in d = 3 di- 
mensions while Ul ^ 0.85Q in d = 2 dimensions 
[111 ]). Consequently, plotting f/4 versus rj^ for 
different L one expects a family of curves that 
intersect at 1]^ = r^p „ in a common intersection 
point, if L is large enough so that corrections to 
finite size scaling are negligible, and using this 
method (or an analogous reasoning [112j for the 
moment ratio M = (m^)/(|m|)^ [8OI |8T] which 
should yield a universal intersection in c? = 3 at 
M* = 1.239 [113]) one finds the estimate of r]^^^^ 
included in Figs. [H HI 

A further consequence of finite size scaling [28l 
[291 [301 [93 EH [991 [M [M [103] is the fact that 
the moments {\m\'') are homogeneous functions 
of the two variables L and t = rf^/rf^f,v ~ 1; 
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'') = L-^PI^M 
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where /3 and u are the critical exponents of the 
order parameter Mc and correlation length ^, 
respectively. 



M, = Bt^, i = ir 



(9) 



In Eq. (jHj), M.k{,x) is a scaling function, and B 
and ^ are critical amplitudes. For the univer- 
sality class of the d = 3 Ising model [114] . the 
exponents are [23, IllSj 1116] 



/3 ^ 0.326, 



z/ 



0.630, 



(10) 



which differ from the corresponding mean-field 
results [114] 



/3mf = 1/2, z/MF = 1/2 



(11) 



Taking the estimates for the exponents 
[Eq. (TTOj) ] and using rj^^^j. = 0.765 we can re- 
plot the data of Fig. H] in scaled form (Fig. [5]), 
and indeed the data collapse rather well on a 
master curve, as implied by Eq. ([H]). If we use 
Eq. (TTT]) instead, no such data collapsing is ob- 
tained. This result shows that finite size scal- 
ing holds, and the AO model also falls in the 
d = 3 Ising universality class, as one might have 
expected. Moreover, the straight line behavior 
seen on the log-log plot for large x not only im- 
plies that the data indeed are compatible with 
the power law, Mc = Bt^, but also the critical 
amplitude B can be estimated with reasonable 
accuracy, B = 0.27 ± 0.02 [53! • This power law 
actually has been included in Fig. H] for 1]^ near 
r^p^^j,. It results from Eq. (jS]) as the asymptotic 
behavior for L —>■ oo. 

As is evident from the insert of Fig. [H the fluc- 
tuations that are ignored by mean-field theory 
j75] have two effects: one effect is that the criti- 
cal point T^p^j. is shifted upward (the compatibil- 
ity of the colloid-polymer mixture is enhanced), 
and the coexistence curve is flattened near the 
critical point [according to mean-field theory, 
Eq. fITT]) . it is a simple quadratic parabola]. 

A similar discussion can be given for the inter- 
facial tension, 7lv (Fig. [3]), which is found to 
vary as [83] 



^Lv = 7t'', /i = 1.26, 7^0.26±0.02, (12) 



while mean-field theory would imply /x = 3/2 
[114] . Vink et al. [83] have also analyzed the 
critical behavior of susceptibilities at both sides 
of the transition and studied the rectilinear di- 
ameter ?7d, as well as a few critical amplitude 
ratios. All these analyses did confirm the Ising 
character of the transition, indicating that the 
Ising critical region in fact is remarkably wide. 
Mean-field theory [75] is only reliable very far 
away from criticality. 



3 Confinement by Sym- 
metric Walls: Evi- 
dence for Capillary- 
Condensation-Like Be- 
havior 

In this section we consider colloid-polymer mix- 
tures in a L X L X D geometry, where confine- 
ment is effected by two identical walls a distance 
D apart. In the simulations, we apply periodic 
boundary condition in the x and y-directions 
parallel to the walls, and again the strategy will 
be to carry out an extrapolation to the thermo- 
dynamic limit via a finite size scaling analysis. 
If one simply uses hard walls for both colloids 
and polymers, as done in [61j, one encoun- 
ters a very pronounced depletion attraction be- 
tween the colloids and the walls, giving rise to a 
very strong "capillary condensation" -like shift 
[TU| [TT] of the coexistence chemical potential 
/icoex(^p) of the colloids. It is hence convenient 
to apply in addition a square-well repulsive po- 
tential 

U^^{h) =e, R^<h< 2i?c, U^^h > 2i?c) = 0, 

(13) 
with h the distance of a colloidal particle from 
the (closest) wall. Of course, Uc^{h < R^) = cxo 
and Up^{h < Rp) = cxo, since neither colloids 
nor polymers are allowed to penetrate into the 
wall. 

If one considers very large e, colloids are ex- 
cluded from the close vicinity of the walls, and 
an effective attraction of the polymers to the 
walls would result. As a consequence, "capillary 
evaporation" is expected rather than "capillary 
condensation" (i.e., close to phase coexistence 
in the bulk the capillary prefers the vapor-like 
phase rather than the liquid-like phase of the 
colloid-polymer mixture). Schmidt et al. [SH] 
presented a (somewhat qualitative) evidence for 
this phenomenon. 

Since the finite size thickness D limits growth 
of the correlation length ^ of volume fraction 
fluctuations near the critical point in the z- 
direction perpendicular to the confining wall, a 
divergence of ^ as described in Eq. 1^ is only 



possible along the x- and y-directions parallel to 
the walls. Therefore, the phase transition which 
can take place is a phase separation in lateral di- 
rections (x, y) only, between colloid-rich and col- 
loid pure phases. As a consequence, ultimately 
this transition should belong to the universality 
class of the two-dimensional Ising model [114] . 
and the critical exponents are 



P = l/t 



H = l, 



(14) 



instead of those quoted in Eqs. fITU]) . flT^ . How- 
ever, this two-dimensional critical behavior pre- 
vails only when ^ has grown to a size much larger 
than D: if ^ -C Z) the behavior is still close to 
three-dimensional, and when ^ and D are of the 
same order a gradual crossover between the two 
types of critical behavior occurs. 
These crossover phenomena make the analysis 
of the simulations somewhat more difficult. For 
any finite value of L the "raw data" estimates 
for r]^ , 7]^ are qualitatively similar, irrespective 
of D (Fig. [H]). Again pronounced "finite size 
tails" occur for these estimates in the vicinity 
of r/p^j,, i.e., for any finite L one finds that rj^ 
and T]^ as estimated from the peak positions of 
P(?7c) fail to merge at r]^^^,., but rather continue 
further into the one-phase region, as in the bulk 
(Fig. H]). When one then plots f/4 vs. r]^ for 
different choices of L, searching for a universal 
intersection point, one rather finds that the in- 
tersection points are somewhat scattered over a 
region of values for 1]^ (Fig. [7]). In addition, this 
intersection does occur neither at the theoreti- 
cal value for U* for the d = 2 universality class 
nor at the U* for d = 3, but rather somewhere 
in between. These findings are a consequence 
of the gradual crossover in critical behavior al- 
luded to above. While for D = 3 both u and 
U* are rather close to the theoretical d = 2 
values, for D = 10 both z/ and U* are about 
half way between the d = 2 and d = 3 values. 
However, these numerical results do not have 
any fundamental significance; they only mean 
that the larger D the closer rj^^^j. needs to be ap- 
proached, to be in the region where ultimately 
^ ^ D and hence the correct asymptotic criti- 
cal behavior (which is always two-dimensional, 
for any finite value of D) can be seen. 
For the case D = 5, e = a. very careful 



analysis has been performed [61], applying a 
novel variant of finite scaling which does not 
imply any bias on the type of critical exponents 
[TTTl fnsl fTT9] . Fig. E shows that the resulting 
order parameter can be fitted over some range 
indeed by an effective exponent PeS = 0.17, 
which is in between the d = 2 and d = 3 values 
(0.125 < /5eff < 0.326), but a correct interpre- 
tation of this finding is that a log-log plot of 
the order parameter vs. t exhibits a slight cur- 
vature, spread out over several decades. Only 
for t -^ can the d = 2 value (/3 = 0.125) be 
expected to be seen; for larger t the slope PeS on 
the log-log plot increases systematically (but it 
does not reach the d = 3 value, since for t > 0.1 
noncritical saturation effects come into play). 
While the critical behavior of thin films is con- 
trolled by the d = 2 critical exponents, a differ- 
ent answer results when one considers the shift 
of the critical point relative to the bulk [33t[34]: 
this shift is controlled by three-dimensional ex- 
ponents only, namely 



oc D-^/^ D^oo (15) 

OC ^(^1-^)/'^, D^oo. (16) 



Here A ^ 1.56 [971 [ffl [IIE] is the so-called 
"gap exponent" which characterizes the bulk 
equation of state near criticality, and Ai ^ 0.47 
[TM fT2n [T221 [T23] its surface analog. Figure 
|9] shows that the AO model is compatible with 
these predictions, even though only rather small 
film thicknesses were accessible to the simula- 
tion (the largest film thickness included in Fig. [9] 
is only for D = 10 colloid diameters). 
Note that for r]^ > r]^,^^ the asymptotic behavior 
of the shift of the colloid chemical potential at 
phase coexistence is not given by Eq. (IT6l) . but 
by the simpler "Kelvin equation" [TOl 135] 
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Fig. [lO] shows that the data of Vink et al. [62] 
are compatible with this equation as expected. 
We emphasize that in mean field theory one 
could not discuss the crossover between two- 



and three-dimensional critical behavior, since 
/9mf = 1 irrespective of dimensionality [114] , 
and also the mean-field predictions for the shift 
of the critical point [Eqs. (fTSll . (fT6ll ] would be 
different from what was observed [S2], since 
l/py^^ = 2 instead of l/u ^ 1.59, and (A^^^ - 
A^^)/z/MF = -2 instead of (Ai - A)/z/ ^ 
— 1.73. However, mean-field theory does repro- 
duce the Kelvin equation, Eq. flT7|) . and in any 
case mean-field results for our confined films 
would be desirable. We note that some mean- 
field results as well as Monte Carlo results are 
available for g = 1 [581 EFlj; however, the accu- 
racy of these Monte Carlo data was too limited 
to allow for a comprehensive test of theoretical 
predictions, as reviewed above, and hence these 
studies EHl ESI are not discussed further here. 



We conclude this section by discussing the 
structure of the coexisting phases in the thin 
film in more detail. Already snapshot pictures 
[Fig. [11] show that in the ^-direction the com- 
position can be inhomogeneous. For e = 
colloidal particles are enriched at the walls in 
the vapor- like phase, while for e = 2 poly- 
mers are enriched at the walls in the liquid- 
like phase. However, it would be wrong to con- 
sider these enrichment layers as wetting layers 
[211 [251 [261 [27] : wetting layers are macroscopi- 
cally thick, and cannot occur in a thin film ge- 
ometry [TT]. One should also recall that wet- 
ting at the surface of a semi-infinite system oc- 
curs at bulk coexistence, while coexistence in 
the thin film deviates from bulk coexistence {cf. 
Eqs. (dnD, ([HD, and Figs. [Si, [ini}. Fig. [H shows 
density profiles across the thin film for several 
typical choices of parameters. One recognizes 
that the colloid density in the liquid -like phase 
near the walls shows a pronounced layering ef- 
fect, while the polymer density in the vapor-like 
phase lacks a corresponding effect. This finding 
is expected, since layering is a consequence of 
the repulsive interactions among the particles. 
While for D = 10 the films do reach homoge- 
neous bulk-like states in their center, for Z^ = 3 
and D = 5 (not shown) the behavior stays in- 
homogeneous throughout the film. 



Confinement by Asym- 
metric Walls: Evidence 
for an Interface Localiza- 
tion Transition 



By asymmetric walls one can realize a situa- 
tion that one wall attracts predominantly col- 
loids and the other wall attracts polymers. As 
discussed in the previous section, hard walls at- 
tract colloids via a depletion mechanism; but 
coating the wall by a polymer brush under 
semidilute conditions, one may cancel this de- 
pletion attraction partially or completely, and 
also reach a situation where polymers get at- 
tracted to the wall. This situation is qualita- 
tively modelled by a step potential of height 
e, acting on the colloids only {Eq. (TT3|) }. An 
asymmetric situation occurs e.g. if the left wall 
is a hard wall but on the right wall the ad- 
ditional potential described by Eq. (1131) acts, 
see Fig. ( 1T5]) . In drawing the schematic phase 
diagrams, we have assumed that for a semi- 
infinite system colloid-polymer mixtures exhibit 
complete wetting [TH |2ll |25l |26l |27| over a 
wide range near the critical point, namely for 
Vl,cv < Vl < Vl,^, while for r]^ > r]^^^ "in- 
complete wetting" (i.e., a nonzero contact an- 
gle of a droplet) would occur. This assump- 
tion is corroborated by density functional cal- 
culations jl2^ and Monte Carlo simulations 
|1U51 1125] . Since no "prewetting transition" 
[m [21 [221 [251 [27] was found, the wetting transi- 
tion presumably is of second order, and this was 
assumed drawing the phase diagrams of Fig. [131 
since this greatly simplifies the theoretical anal- 
ysis. For the case of symmetrical mixtures of 
long fiexible macromolecules, the infiuence of 
prewetting phenomena on the phase diagram of 
thin confined films has been thoroughly investi- 
gated [IIl[39l[65l[66l[671[68l[69l[70l[ni[I2], and 
it has been shown that typically a phase dia- 
gram with two critical points and a triple point 
can be expected. 

For asymmetric walls an interface localization 
transition may occur, and this situation is ex- 
plained qualitatively in the right part of Fig. [121 
If the strength of the attraction of the colloids 
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of the left wall is of the same order as the 
strength of the attraction of the polymers to the 
right wall, 1]^^^^ and ^p'^^* will be rather close 
to each other and both exceed rj^^^j. distinctly. 
Then for r^p > t]'^ j. the left wall will always 
be coated with colloids, the right wall will al- 
ways be coated with polymers. In other words, 
we expect an interface between the colloid-rich 
phase on the left and the polymer-rich phase 
on the right. When fi is small enough {i.e, 
/Wcocx(oo)— /i is large enough} most of the system 
is in the polymer-rich phase (shown schemati- 
cally as Bllb in Fig. [T3|) but when /z increases 
a transition takes place to a state where most 
of the film is in the colloid-rich phase (state 
Blla). For r]^ > rf^cvi^) ^^is transition is 
a sharp (first-order) phase transition, i.e. the 
interface jumps from a state localized at the 
left wall to a state localized at the right wall. 
For if = ri^^j.{D) this transition is of first or- 
der, while for r^Jj < rf ^^{D) the transition is a 
smooth gradual transition (near the broken line 
in Fig. [T3l) . Note, however, that this transition 
becomes sharper and sharper as D increases, 
but a true phase transition appears only in a 
discontinuous manner in the limit D ^ oo [TT]: 
then the broken line in Fig. [T3] coincides with the 
hue yu = /icoex(oo) cudiug at r^p^.^ , and rf^^^AD) 
does not converge to rf„^^ but rather we have 

in Fig.[I3D. 

For the states along the broken curve in Fig. [T3] 
the system is essentially inhomogeneous, there 
exists a thick domain of colloid-rich phase in 
the left part of the film, and a thick domain 
of polymer- rich phase in the right part, sepa- 
rated by a "delocalized" interface in the cen- 
ter of the film [IH EHl EZl ESI EH]. Snapshot 
pictures of the system indeed readily confirm 
such a scenario (Fig. [H]), as well as the density 
profiles across the thin film (Fig. [T5|) . Fig. [T5b 
shows the profile for rf^ = 0.7 < rj^cr^ i-6- a 
state in the one phase region of the bulk. One 
recognizes that the colloid concentration is en- 
hanced near the hard wall, as expected from 
the depletion attraction. Near the other wall at 
z = D = 10, the colloid concentration is some- 
what depressed, but the polymer concentration 
is clearly enhanced. But in the center of the thin 



/p,cr 

^p,w (fo'^ the situation drawn 



film both profiles are roughly constant, as ex- 
pected for bulk-like behavior. In fact, for large 
D we expect that the surface enhancement (or 
reduction, respectively) decays with z according 
to an exponential relation, exp(— /i/^), where h 
is the distance from the closest wall and ^ the 
bulk correlation length |120j . Figure [T5b shows 
the profiles at r]^ = 0.95, which exceeds the bulk 
critical value ^pcr? but still is smaller than the 
critical value rjp^cri^) of ^^e confined system. 
Now, the profiles are very different from those 
of Fig. [T5h : phase separation in a colloid-rich 
and a polymer-rich phase has occurred, with 
the interface position (estimated from the in- 
flection point of the polymer volume fraction 
profile rip{z), for instance) being located in the 
center of the film. This situation corresponds to 
the snapshot in Fig. [TDo. The interfacial profile 
resembles that of an interface between bulk co- 
existing phases, broadened by capillary waves 
PI Eg. Finally, the cases ril = 1.2 (Fig. [I5t,d) 
refer to the two-phase region of the film. The in- 
terface either is located near the hard wall, cor- 
responding to the polymer-rich phase of the film 
(Fig. [T5b : this case corresponds to the snapshot 
shown in Fig.fT^). or near the wall that attracts 
the polymers (Fig. [T^). Note that along the 
transition line drawn schematically in Fig. [T3] 
(right part), there occurs lateral phase separa- 
tion between the states corresponding to these 
two types of profiles. Fig. \Wc and Fig. \Wd . 
which hence can coexist with each other in a 
thin film (and then are separated by an inter- 
face running from the right wall towards the left 
wall). 

Figure [16] shows the corresponding phase di- 
agrams for two film thicknesses, D = 5 and 
D = 10, varying the strength e of the poten- 
tial [Eq. (TT^ ] at the right wall. One sees that 
with increasing e the critical points and the 
whole coexistence curves are shifted upwards, 
to rather large values of rj^. This shift is con- 
sistent with the qualitative phase diagram of 
Fig. dni (right part). Of course, one must again 
recall that in Fig. [16] we show "raw Monte Carlo 
data" for one choice of L only, and hence pro- 
nounced finite size tails near rjpcriD) are appar- 
ent, as discussed for the case of capillary con- 
densation already (Fig. [Sj). The critical points 
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were again estimated from the cumulant inter- 
section method. Although strong corrections to 
scahng are present, the conclusion can be drawn 
[63j that the critical behavior of the interface lo- 
calization belongs to the d = 2 Ising universality 
class, as expected [3Sl EH EH] • 
Figure [17] shows the phase diagram for D = 5 
and various choices of e in the grand-canonical 
representation (this is the counterpart of Fig. [9ti 
for capillary condensation, where D was varied 
for e = 0, while here we study the variation with 
e at fixed D). One can see that by increasing 
e the coexistence curves of the thin film move 
closer towards the bulk coexistence curves, and 
for e = 2.5 the deviation from the bulk indeed 
is very small, but the critical point is strongly 
shifted (from rj^^^j. = 0.766 in the bulk to ?7p „ = 
1.106 in the thin film). 

This behavior is qualitatively similar to what 
has been found for the Ising ferromagnet with 
competing surface magnetic fields [HI EH] , the 
generic model for which the interface localiza- 
tion transition was studied for the first time. 
Note that the curve for e = in Fig. [T7] 
represents capillary condensation (and a simi- 
lar conclusion applies to the case e = 0.5 as 
well). Figure [T7| implies that varying e one 
can completely smoothly cross over from cap- 
illary condensation- like behavior to interface 
localization-like behavior, when e is increased. 
In view of the qualitative description of Fig. [12] 
this is somewhat surprising: in the capillary 
condensation transition, the two liquid-vapor 
interfaces bound to the walls annihilate each 
other, the slit pore gets almost uniformly filled 
with liquid. In the interface localization transi- 
tion, one has an interface on both sides of the 
transition, it just has jumped at the transition 
from one wall to the other. 
How can one then reconcile Figs. [13] and [17] 
with each other? The clue to the problem is, of 
course, that the picture of the states in Fig. [13] 
is far too simplified, it ignores the variations 
of the densities close to the wall. Therefore 
the states with "interfaces bound to walls" are 
a simplification, which lose its meaning when 
the "phase" in between the interface and the 
wall can no longer be clearly identified with 
bulk-like properties (as is actually the case, see 



Fig. [TB]) . While one can clearly imagine to trans- 
form the left phase diagram of Fig. [13] into the 
right one by smooth changes, one should not 
take the sketches that illustrate the character 
of the phases too literally. The failure of these 
sketches, however, also means that one must be 
careful with all approaches where wetting phe- 
nomena and interface localization transitions 
[21 1251 [261 [23 [361 [an [38] are simply described 
in terms of the "interface hamiltonian" picture, 
since according to this description the distance 
i of the interface from the wall is the single de- 
gree of freedom (on a mean-field level) left to 
analyze the problem. 



5 Summary and Outlook 

In this brief review we have emphasized that 
confinement has very interesting effects on soft 
matter systems, both with respect to the struc- 
ture and the phase behavior of these systems. 
Of course, confinement also has very interest- 
ing consequences on the dynamics of soft mat- 
ter systems (see e.g. [El [El [151 [ISl [IH [IHl [H] 
for recent discussions) , but this aspect has been 
completely outside of the focus of our review. 
We also have focused on the case which we con- 
sider to be the simplest case, confinement be- 
tween two fiat and ideally parallel walls a fi- 
nite distance D apart. Practically more impor- 
tant, of course, is the confinement in random 
porous media [IDl [HI [201 [211 [221 [23] • However, 
in this case the random irregularity of the con- 
fining geometry is a serious obstacle for a de- 
tailed understanding. There is ample evidence 
(both from experiment [1261 I127J and simula- 
tions [I2Hl[I2Sl[ISn]) that the liquid-vapor type 
phase separation or demixing of binary fiuid 
mixtures under such confinement is seriously 
modified, but the character of this modifica- 
tion has been under discussion since a long time 
[I2S1 [123 [I2HI [I2S1 11201 [Ml [I32]. De Gennes 
|131j argued that due to the random arrange- 
ment of the pore walls (which prefer one of the 
coexisting phases over the other) the problem 
can be mapped to the random field Ising model 
|1331I134] . While for a long time the existing ev- 
idence [I2S1 [123 [I2H1 112211120] was inconsistent 
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with this suggestion, in recent work on colloid- 
polymer mixtures confined by a fraction of col- 
loids that are frozen in their positions and do 
not take part in the phase separation, evidence 
for the random field Ising behavior was obtained 
[132] . Note that then it is necessary to reach a 
regime where the correlation length has grown 
to a large enough distance, much larger than the 
characteristic linear dimension of the confining 
particles. 

Another case, that has received ample consider- 
ation in the literature (see [12] for further ref- 
erences) but was disregarded here, is the con- 
finement in a quasi-one-dimensional cylindri- 
cal geometry. In this case again the structure 
is typically inhomogeneous in the radial direc- 
tion perpendicular to the walls of the cylin- 
der. The correlation length can grow indefi- 
nitely only in one direction, along the cylinder 
axis, however. Therefore the phase transition 
to this laterally segregated state is a gradual 
(rounded) transition only, and even for condi- 
tions where in the bulk the system is strongly 
segregated (with an interfacial tension 7 be- 
tween the coexisting phases which is not small in 
comparison with ksT), there is no macroscopic 
phase separation possible: rather one predicts 
that only phase separation into domains of fi- 
nite size can occur (cross sectional area A of the 
cylinder and length id of the domains), where 
|135[ I136J In id oc A'y/kBT. This happens be- 
cause thermal fiuctuations prevent the estab- 
lishment of true long range order by the sponta- 
neous generation of transverse interfaces (across 
the cylinder). This has the consequence that 
also, in principle, the capillary condensation or 
evaporation transitions in cylindrical geometry 
are, in full thermal equilibrium, not perfectly 
sharp but rounded. In practice, this effect often 
is masked by nonequilibrium phenomena (pro- 
nounced hysteresis occurs!) and hence we are 
not aware of careful studies where this round- 
ing has been demonstrated. Note that mean- 
field treatments [T2] miss such fiuctuations ef- 
fects, of course. Also, when one considers cylin- 
drical geometries with '"competing walls" (e.g. 
a cylinder with a square cross section, where 
the upper walls prefer one phase and the lower 
walls prefer the other phase |137[ I138J ) the in- 



terface localization transition in such a "double 
wedge" -geometry is rounded and not sharp, as 
simulations show where one considers the limit 
that the length of the cylinder gets macroscopic 
while its cross section stays finite |138j . Since 
recently it has become possible to create artifi- 
cial cylindrical nanochannels with diameters be- 
tween 35 and 150 nm |139j, it would be interest- 
ing to study phase separation in such nanochan- 
nels experimentally as well. 
Using the AO model of colloid-polymer mix- 
tures as an example that is well suited for sim- 
ulation studies [nil E21 E3], we have discussed 
simulation evidence for the theoretical concepts 
on capillary condensation and interface localiza- 
tion transitions [IIl[33l[32[35l[36l[371[38]. In 
particular, the predictions for the shift of the 
critical point have been found to be compati- 
ble with the simulation results, and it was also 
argued that the critical behavior of the lateral 
phase separation in the thin film has the char- 
acter of the two-dimensional Ising model (al- 
though in practice one is mostly in a crossover 
region where "effective exponents" in between 
the d = 2 and d = 3 limits apply, which do not 
have a deep theoretical significance). Clearly, 
it would be nice to have also experiments that 
confirm the findings of theory and simulation on 
the phase behavior of confined fiuid mixtures. 
One crucial assumption of the work reviewed 
here was that the wetting transitions (that oc- 
cur in the limit when the film thickness D di- 
verges to infinity) are of second order [11]. If 
first-order wetting occurs, much more compli- 
cated phase diagrams under confinement result 
|11] . To some extent, this problem has been 
worked out for symmetrical polymer blends un- 
der confinement [65l[66l[67l[68li9l[70l[ni[I2], 
and we refer the reader to these papers for de- 
tails. In particular, it also should be possible 
to realize situations in between capillary con- 
densation and interface localization transitions 



Finally, we draw attention to more compli- 
cated ordering phenomena under confinement. 
A problem that has found a lot of attention is 
the effect of confinement in thin films on the 
blockcopolymer mesophase ordering |14Ul I14H 
I142[ 1143] . For more or less symmetric compo- 
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sition of a diblock copolymer, the mesophase erker, H. Lowen, M. Miiller, M. Schmidt and P. 
observed in the bulk is a lamellar ordering [5^ . Virnau for many stimulating discussions. 
The question that is then discussed in the lit- 
erature (both experimentally and theoretically, 
see |14Ul 11411 1142] for further references), is 
whether the lamellae are oriented parallel or 
perpendicular to the confining walls, and tran- 
sitions in the number of lamellae that fit into 
the thin film, etc. (remember that the lamella 
thickness depends on temperature, chain length 
of the polymers, and other control parameters 
[57]). For asymmetric compositions AfBi^j of 
a diblock copolymer, however, already in the 
bulk melts other mesophases appear, such as 
hexagonal patterns of A-rich cylinders in B-rich 
background, or cubic structures, where A-rich 
cores of micelles form a periodic lattice in the 
B-rich background, or vice versa |57|. For tri- 
block copolymers, many much more complex 
mesophases occur, and the question how all this 
self-assembly of block copolymers is affected by 
confinement due to walls is still under both 
theoretical and experimental investigation (see 
e.g. tl43j). 

Other very interesting confinement effects in 
soft matter occur when orientational order is 
involved, e.g. when a colloidal dispersion un- 
dergoing a transition from isotropic to nematic 
phases is confined by walls (see e.g. [144[ 11451 
11461 1147] ). Confinement may enhance the ne- 
matic ordering tendency ("capillary nematiza- 
tion" [144] is the analog phenomenon of capil- 
lary condensation), but one needs also to take 
into account the tensor character of the order 
parameter of liquid crystals. Thus near a wall 
a biaxial character of the ordering occurs even 
when in the bulk the ordering is uniaxial. Also 
the boundary conditions at the walls can be en- 
visaged such that one wall prefers parallel and 
the other wall perpendicular alignment, lead- 
ing to a tilted structure of the ordering across 
the film [145] . These remarks are by no means 
intended as an exhaustive discussion, but just 
want to draw the attention of the reader to this 
wealth of interesting problems. 
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Figure 1: Phase diagram of the AO model with 
a size ratio q = Rp/Rc = 0.8 plotted in the 
('7p,''?c) plane, showing the coexistence curve 
for phase separation into a polymer-rich phase 
(left) and a colloid-rich phase (right), according 
to Monte Carlo (circles) and the free volume 
mean-field theory of Lekkerkerker et al. (full 
line) [75]. The Monte Carlo data were obtained 
using a simulation box with linear dimensions 
Lx = Ly = 16.7, Lz = 33.4 (some data for a 
smaller box with Lx = Ly = 13.3, Lz = 26.5 
are also included as crosses). The solid square 
shows the estimate for the critical point as ob- 
tained from a finite size scaling analysis of the 
Monte Carlo data, cf. text. The insert shows the 
same data, in the so-called "reservoir represen- 
tation" where ?7p rather than rip is used as a vari- 
able. The one-phase region where colloids and 
polymers are fully miscible is always shown in 
the lower parts of these diagrams, below the co- 
existence curves. From Vink and Horbach 
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Figure 2: Logarithm of the probability P{rjc) 
of observing a colloid packing fraction rjc for an 
AO mixture with q = 0.8 at coexistence for sev- 
eral values of ?7p as indicated. The simulations 
were performed in a box with linear dimensions 
Lx = Ly = 16.7 and L^ = 33.4 using periodic 
boundary conditions. Note that the distribu- 
tions are not normalized, and that for rj near 
Vd = (^7^ + ^^)/2 the distribution is essentially 
flat, almost independent from rjc, so that a free 
energy barrier AF is well-defined. From Vink 
and Horbach IST]. 
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Figure 3: Reduced interfacial tension 7* = 
AR^'-^YM/ksT for an AO mixture with q = 0.8 
as a function of the difference in packing frac- 
tion between the coexisting phase, rj^ — r]^. The 
inset shows 7* vs. ?7p. Open circles are the sim- 
ulation results, while the full curves are density 
functional theory predictions. From Vink and 
Horbach EH 
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Figure 4: Close-up of the phase diagram of the 
AO model for q = 0.8 in the critical region, 
showing only the range 0.75 < r]^ <= 0.80 for 
the polymer reservoir packing fraction. Results 
for four different choices of the linear dimension 
L oi L X L X L simulation boxes are shown as 
broken and dotted lines. The extrapolation to- 
wards the thermodynamic limit, L = 00, also is 
included, and the estimate for the critical point 
'7pcr — 0.766 ± 0.002 is shown by the error bar. 
From Vink et al. I82] 
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Figure 5: Scaling plot for the order parameter 
Mc = {\Tn\) of the colloid-polymer mixture in 
the bulk, for the AO model with q = 0.8, using 
only data for t > 0. The quantity L^^'^ Mc is 
plotted vs. tL^^'^, choosing logarithmic scales. 
The straight line is a power law y = 0.27a:^, 
with f] = 0.326. From Vink et al. 




Figure 6: Coexistence curves of the AO model 
with q = 0.8 in the bulk (full curve without 
data points, representing an extrapolation to 
the thermodynamic limit) and for thin films 
confined by symmetric walls, choosing Eq. (fT3ll 
with e = 0.5 as a wall potential. Open circles 
denote data for D = 3, L = 18 (all lengths 
here are measured in units of the colloid diam- 
eter, (Tc = 2Rc). Open squares denote data for 
D = 5, L = 20; open diamonds are for D = 7.5, 
L = 30 , and open triangles for D = 10, L = 30. 
The broken curves denote the corresponding co- 
existence diameters, rjd = (r/^ + ^^)/2 and the 
full symbols denote estimates of the correspond- 
ing critical points for the various choices of D. 
From Vink et al. i62l. 



15 




0.885 0.890 0.895 0.900 



Figure 7: Moment ratio [4 for a film of thick- 
ness D = 5 plotted versus the polymer reser- 
voir packing fraction ?7p and several choices of 
the lateral linear dimension, L = 12.5, 15, 17.5, 
20, 22.5, and 25. The intersection points al- 
low to estimate the critical point as 7]^,^^. = 
0.892 ± 0.002. The inset shows a log-log plot of 
the slope Yi = dU^/dr]^ at ?7p^j. versus L, test- 
ing the prediction that Yi oc L^^'^. The straight 
line in the inset corresponds to an (effective) 
exponent Ugg = 0.933. The broken horizontal 
straight lines in the main part indicates the the- 
oretical values for the universal value U* of ^4 
at the intersection point, for both the d = 2 and 
the d = 3 Ising model, respectively. From Vink 
et al. [62] . 
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Figure 8: Order parameter A(= Mc) of the con- 
fined AO model with q = 0.8 plotted vs. rj^, 
choosing D = 5 and e = 0. The curve through 
the simulation data in the main frame is a fit to 
A = fieffi^"«, choosing r/^„ = 0.9223. The inset 
shows the same data as a function of the relative 
distance from the critical point, t = rf^/rf^^cv ~ !> 
on double logarithmic scales; broken straight 
lines illustrate the two-dimensional and three 
dimensional Ising exponents. From Vink et 
al. EH. 
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Figure 9: a) Coexistence curves of the AO 
model with q = 0.8 in the grand-canonical rep- 
resentation where the chemical potential n'^°^^ 
of the colloids at the coexistence curve is plot- 
ted vs. the polymer reservoir packing fraction 
7]p. The bulk result [D = cxd) is shown as a 
full curve, while the broken curves show the co- 
existence curves for confined films for several 
thicknesses D, as indicated. The symbols mark 
the corresponding critical points, b) Shift of 
the critical coexistence colloid chemical poten- 
tial plotted vs. £)(Ai-A)/!. Equation (^ im- 
plies a straight line, as indicated by the dashed 
curve, c) Shift of the critical polymer reservoir 
packing fraction, Ar]p,,j.{D), plotted vs. D'^/" . 
The dashed lines indicate that Eq. flTSl) holds. 
From Vink et al. 
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Figure 10: Test of the Kelvin equation. 
The chemical potential difference A/i'^°''^(D) 
{Eq. (ITTj) } is plotted vs. D~^ for three val- 
ues of ?7p, chosen well above the critical values 
rf ^^{D\ Broken straight lines show that the 
data are compatible with the Kelvin equation. 
From Vink et al. [62]. 



a)e = 0.0, Ti^=l.l 



b) e = 2.0, Ti^=l.l 



" •-; r , si.'* . z 

vapor 






liquid 




vapor 




liquid 



Figure 11: Snapshot pictures of coexisting 
phases for the colloid-polymer mixture with q = 
0.8, D = 10, 7]l = 1.1, for e = (a) and e = 2.0 
(b) Colloidal particles are in green, polymers in 
blue (the size of the polymers is rescaled to al- 
low a clearer view). 
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Figure 12: Colloid density profiles obtained 
in thin films at rj^ = 1.1, for two values of 
the film thickness D, and several values of the 
colloid-wall parameter e as indicated. Frames 
a) and b) show profiles obtained for D = 10, 
on the vapor and liquid branch of the coex- Virgiliis et al. 
istence curve, respectively. Frames c) and d) 
show the corresponding profiles for thickness 
D = 3. Note the jumps in the colloid density at 
z = 0.5, 1.0, 2.0, 2.5 caused through the jumps 
of the potential at z = Re and z = 2Rc, respec- 
tively. The insets represent density profiles of 
the polymers. From Vink et al. |62j . 



Figure 13: Schematic phase diagram of a 
colloid-polymer mixture confined between two 
parallel walls a distance D apart, in the grand- 
canonical ensemble where the polymer reservoir 
packing fraction i]^ is used as ordinate and the 
difference between the chemical potential of the 
colloids at bulk phase coexistence /icoex(-D = oo) 
and the actual colloid chemical potential is used 
as abscissa. Thus, phase coexistence in the 
bulk occurs along a vertical straight line at 
Aicoex(oo) — fi = 0. The left part shows the case 
of symmetric walls, the right part asymmetric 
walls. States Al and All coexist along the cap- 
illary condensation transition line, states Blla 
and Bllb coexist along the interface localiza- 
tion line, while state BI exists along the broken 
curve (which represents a line of rounded transi- 
tions). In the limit D ^ oo, which corresponds 
to an infinite system but bounded by walls both 
on the left and the right side, wetting transi- 
tions occur, that are rounded off for finite D. 
In the symmetric situation, the wetting transi- 



tions of both walls coincide at ?7p,„, while in the 



asymmetric situation r/''''''§^* 



/p,w 



/p,W5 

_! ^r,left 
/ 'P,w ■ 



From De 



a)e = 2.5, D = 10,Ti =0.18, Ti' = 0.7 : 



a) 11;;= 1.1 



£0= 2.0 




b) 11;;= 0.9 
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£0= 2.0 




£0 = 0.0 

Figure 14: a) Snapshot picture of the polymer- 
rich phase, for D = 10, q = 0.8, rj^ = 1.1. The 
lower wall (at z = 0) has a wall potential param- 
eter energy £0 = and hence attracts colloidal 
particles (shown in green), while the upper wall 
(at z = D) has a wall potential energy parame- 
ter En = 2, attracting polymers (shown in blue), 
b) Snapshot picture of the same system as in a), 
but for ?7p = 0.9 showing a state with a delocal- 
ized interface. 
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Figure 15: Colloid concentration profiles ric{z) 
and polymer concentration profiles rip{z) as a 
function of z for a thin film with asymmetric 
walls (hard wall at 2; = 0, while for the other 
wall at z = D = 10 the potential Ucw{h) acts see 
Eq. (TT3!) . with e = 2.5). Profiles were obtained 
at ri, = 0.18, ril = 0.70 (a) , r], = 0.18, ril = 
0.95 (b), r/c = 0.05, 7]^ = 1.20 (c), and r]^ = 
0.33, r]l = 1.20 (d). For profiles (c) and (d), the 
choices tjc = 0.05, 0.33 roughly correspond to 
the two branches of the coexistence curve, see 
Fig. [16l From De Virgiliis et al. [63] . 
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Figure 16: Coexistence curves for D = 5, L = 
30 (a) and _D = 10, L = 20 (b), using four 
values of e, as indicated. Also the bulk coex- 
istence curve is shown (full curves). Full sym- 
bols mark critical points, the broken lines end- 
ing at these critical points are the coexistence 
diameters. The dotted horizontal straight lines 
mark the values of yf ^^{D). From De Virgiliis 
et al. [631. 




Figure 17: Phase diagram of a thin film of 
colloid-polymer mixtures with asymmetric walls 
in the grand-canonical ensemble, choosing q = 
0.8, D = 5, with ?7p as abscissa and fi as or- 
dinate. Curves show the coexistence potential 
ficoex{D,e) of the colloids. Full curve denotes 
the result for the bulk (note that the bulk crit- 
ical point, ?7pj.j. = 0.766 is off the scale of this 
figure). Full symbols mark the critical points of 
the films. From De Virgiliis et al. ^63] 
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